Evolution of a quantum spin system to its ground state: Role of entanglement and 
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We study the decoherence of two ferro- and antiferromagnetically coupled spins that interact with 
a frustrated spin-bath environment in its ground state. The conditions under which the two-spin 
system relaxes from the initial spin-up - spin-down state towards its ground state are determined. 
It is shown that the two-spin system relaxes to its ground state for narrow ranges of the model 
parameters only. It is demonstrated that the symmetry of the coupling between the two-spin system 
and the environment has an important effect on the relaxation process. In particular, we show that 
if this coupling conserves the magnetization, the two-spin system readily relaxes to its ground state 
whereas a non-conserving coupling prevents the two-spin system from coming close to its ground 
state. 
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I. INTRODUCTION 

The foundations of non-equilibrium statistical mechan- 
ics are still under debate (for a general introduction to 
the problem, see, e.g., Ref. [l|; see also a very recent 
discussion Q and Refs. therein). There is a common be- 
lieve that a generic "central system" that interacts with 
a generic environment evolves into a state described by 
canonical ensemble (in the limit of low temperatures, this 
means the evolution to the ground state). Experience 
shows that this is true but a detailed understanding of 
this process, which is crucial for a rigorous justification of 
statistical physics and thermodynamics, are still lacking. 
In particular, in this context the meaning of "generic" 
is not clear. The key question is how the evolution to 
the equilibrium state depends on the details of the dy- 
namics of the central system itself, on the environment, 
and on the interaction between the central system and 
environment. 

In one of the first applications of computers to a ba- 
sic physics problem Fermi, Pasta, and Ulam attempted 
to simulate the relaxation to thermal equilibrium of a 
system of interacting anharmonic oscillators The 
results obtained appeared to be counterintuitive, as we 
know now, due to complete intcgrability (in the contin- 
uum medium limit) of the model they simulated 0]. 

Bogoliubov Q has considered in a mathematically rig- 
orous way the evolution to thermal equilibrium of a clas- 
sical harmonic oscillator (central system) connected to 
the environment of classical harmonic oscillators which 
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are already thermalized (for a generalization to a nonlin- 
ear Hamiltonian central system with one degree of free- 
dom, see in Ref. 0). Also, for quantum systems this 
"bosonic bath" is the bath of choice, starting with the 
seminal works by Feynman and Vernon[7| and Caldeira 
and Leggett 8] (for a review, see Ref. |9|). On the other 
hand, as we know now, the bosonic environment differs in 
many ways from, say, a spin-bath environment (such as 
nuclear spins) that dominate the decoherence processes 
of magnetic systems at low enough temperatures [To| . 
The evolution of quantum spin systems to the equilib- 
rium state has been investigated in Refs. [ll|, EH O, for 
a very special class of spin Hamiltonians. 

In terms of the modern "decoherence program" quan- 
tum systems interacting with an environment evolve to 
one of robust "pointer states", the superposition of the 
pointer states being, in general, not a pointer state [Til . 
[la ] . The decoherence program is supposed to explain the 
"Schrodinger cat paradox" , that is, the inapplicability of 
the superposition principle to the macroworld. It is con- 
firmed in many ways, indeed, that for the case where 
the interaction with environment is strong in compari- 
son with typical energy differences for the central sys- 
tem classical "Schrodinger cat states" are the pointer 
states. At the same time, some less trivial pointer 
states have been found in computer simulations of quan- 
tum spin systems for some range of the model param- 
eters [TH, 03, EH- In fact, the evolution of quantum 
spin systems to equilibrium is still an open issue (see 
also Refs. Ojjl, Ho|, Recently, the effect of an envi- 

ronment of N 1 spins on the entanglement of the two 
spins of the central system has attracted much atten- 
tion [iEiiiiiii^iiaHiiii^i^isiiiiaS. 

The relationship between the pointer states and the 
eigenstates of the Hamiltonian of central system is of 
special interest for the foundations of quantum statis- 
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tical mechanics: The standard scenario assumes that the 
density matrix of the system at the equilibrium is diago- 
nal in the basis of these eigenstates. Paz and Zurek |3Xj ] 
have conjectured that pointer states are the eigenstates 
of the central system if the interaction of the central sys- 
tem with each degree of freedom of the environment is a 
perturbation, relative to the Hamiltonian of the central 
system. In view of the foregoing, it is important to es- 
tablish the conditions under which this conjecture holds 
and to explore situations in which the interaction with 
environment can no longer be regarded as a perturbation 
with respect to the Hamiltonian of the central system. 

In our Letter [22| , we reported a first collection of re- 
sults for an antiferromagnetic Heisenberg system coupled 
to a variety of different environments. Our primary goal 
was to establish the conditions under which the central 
system relaxes from the initial spin-up - spin-down state 
towards its ground state, that is the maximally entan- 
gled singlet state. We found that environments that 
exhibit some form of frustration, such as spin glasses 
or frustrated antiferromagnets, may be very effective in 
producing a final state with a high degree of entangle- 
ment between the two central spins. We demonstrated 
that the efficiency of the decoherence process decreases 
drastically with the type of environment in the follow- 
ing order: Spin glass and random coupling of all spins to 
the central system; Frustrated antiferromagnet (triangu- 
lar lattice with the nearest-neighbors interactions); Bi- 
partite antiferromagnet (square lattice with the nearest- 
neighbors interactions); One-dimensional ring with the 
nearest- neighbors antiferromagnetic interactions [22[. 

Competing interactions, frustration and glassiness pro- 
vide a very efficient mechanism for decoherence whereas 
the difference between integrable and chaotic systems is 
less important fl8j. Furthermore, we observed that for a 
fixed system size of the environment and in those cases 
for the decoherence is effective, different realizations of 
the random parameters do not significantly change the 
results. However, maximal entanglement in the central 
system was found for a relatively narrow range of the cou- 
plings between the environment spins and the interaction 
between the central spins and those of the environment. 

Having established that the decoherence caused by a 
coupling to a frustrated, spin-glass-like environment can 
be a very effective, it is of interest to study in detail, 
the time evolution of the central system coupled to such 
an environment. In this paper, we consider as a central 
system, two ferro- or antiferromagnetically coupled spins 
that interact with a spin- glass environment. The interac- 
tions between each of the spin components of the latter 
are chosen randomly and uniformly from a specified in- 
terval centered around zero, making it very unlikely that 
there are conserved quantities in this three-component 
spin-glass. For the interaction of the central system with 
each of the spins of the environment we consider two 
cases. 

In the first case, the couplings between the three com- 
ponents are generated using the same procedure as used 



for the environment. In the second case, the central sys- 
tem interacts with the environment via the z-components 
of the spins only. This implies that both the Hamiltoni- 
ans that describe the central system (isotropic Heisen- 
berg model) and the interaction between the central sys- 
tem the environment commutes with the total magneti- 
zation of the central system, hence the latter is conserved 
during the time evolution. In contrast to the naive pic- 
ture in which the presence of conserved quantities re- 
duces the decoherence, we find that the presence of a 
conserved quantity may affect significantly the nature of 
the stationary state to which the central system relaxes. 

II. MODEL 

The model Hamiltonian that we study is defined by 

H = H c + H e + H ce , 
H c = — JSi • S 2 , 

AT-l N 

tf e = -E E E«S } W 

i— 1 j— 2+1 a 
2 N 

^ e = -EEE A S ) w. (i) 

i— 1 j — 1 a 

where the exchange integrals J and Vt^ determine 
the strength of the interaction between spins S„ = 
(Sn,Sv,S*) in the central system (H c ), and the spins 
I n = (1%, 1%, I*) in the environment (i? e ), respectively. 
The exchange integrals A-°^ control the interaction (H ce ) 
of the central system with its environment. In Eq. (fT]), 
the sum over a runs over the x, y and z components 
of spin-1/2 operators S and I. The exchange integral 
J of the central system can be positive or negative, the 
corresponding ground state of the central system being 
ferromagnetic or antiferromagnetic, respectively. 

In the sequel, we will use the term "Heisenberg-like" 
H ce (H e ) to indicate that (0$) are uniform ran- 

dom numbers in the range [-A|J|, A|J|] J|, fi| J|]) 

for all a's and use the expression "Ising-like" H ce (H e ) 

to indicate that A^ y) = (fi^' v) = 0), and that A$ 

( z) 

(fijy) are dichotomic random variables taking the val- 
ues ±A (±fi). The parameters A and determine the 
maximum strength of the interactions. 

The quantum state of central system is completely de- 
termined by its reduced density matrix, the 4x4 matrix 
that is obtained by computing the trace of the full den- 
sity matrix over all but the four states of the central 
system. In our simulation work, the whole system is as- 
sumed to be in a pure state, denoted by \^{t))- Although 
the reduced density matrix contains all the information 
about the central system, it is often convenient to char- 
acterize the state of the central system by other quanti- 
ties such as the correlation functions (5'(t)|Si • S2 1 'I' (£)) , 
(^(t)\SfSI\^(t)), and (^(t)\SfSi\^(t)}, the single-spin 
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TABLE I: The values of the correlation functions (Si • S2), (S'fS'f }, {S1S2), the total magnetization M, the concurrence C and 
the magnetization (Sf) for different states of the central system. 
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magnetizations (*(i)|5f|#(i)), (*(i)|5f |$(t)), and M = 
(#(i)|(Sf + Sf)|*(*)), and the concurrence C{t) 0, 
l34j . The concurrence, which is a convenient measure for 
the entanglement of the spins in the central system, is 
equal to one if the state of central system is unchanged 
under a flip of the two spins, and is zero for an unen- 
tangled pure state such as the spin-up - spin-down state. 
In Table HI we show the values of these quantities for to 
different states of the central system. 

As the energy of central system is given by — (t) | Si • 
S2| X I'(<)), it follows from Table|T]that the four eigenstates 
of the central system H c are given by 



\S) 
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satisfying 



H c \S) = E s \S) , H c |Ti i0 ,_i) = E T |Ti >0 ,_i) 



(2) 



(3) 



where E s = 3 J/4 and E T = - J/4. 

^From Table HI it is clear that the singlet state 15*) is 
most easily distinguished from the others as the central 
system is in the singlet state if and only if (Si • S2) = 
—3/4. To identify other states, we usually need to know 
at least two of the quantities listed in Table Q] For exam- 
ple, to make sure that the system is the triplet state \Tq), 
the values of (Si • S2) and (SfS%) should match with the 
corresponding entries of Table |IJ Likewise, the central 
system will be in the state ||f) if (Si • S2) and M agree 
with the corresponding entries of Table [H 

In general, we monitor the effects of the decoherence by 
plotting the time dependence of the two-spin correlation 
function (Si • S2) and the matrix elements of the density 
matrix. We compute the matrix elements of the density 
matrix in the basis of eigenvectors of the central system 
(see Eq. ©). If necessary to determine the nature of the 
state, we consider all the quantities listed in Table fl] 



The simulation procedure is as follows. First, we se- 
lect a set of model parameters. Next, we compute the 
ground state |</>o) of the environment and, for reference, 
the ground state of the whole system also. The spin-up - 
spin-down state (|TI)) is taken as the initial state of the 
central system. Thus, the initial state of the system reads 
\ty(t = 0))) = |t|) |0o) and is a product state of the state 
of the central system and the ground state of the envi- 
ronment which, in general is a (very complicated) linear 
combination of the 2 N basis states of the environment. 

The time evolution of the whole system is obtained 
by solving the time-dependent Schrodinger equation for 
the many-body wave function \^(t)), describing the cen- 
tral system plus the environment. The numerical method 
that we use is described in Ref . [HJ • It conserves the en- 
ergy of the whole system to machine precision. 

In our model, decoherence is solely due to fact that the 
initial product state 1^(0)) = | Tl) evolves into an entan- 
gled state of the whole system. The interaction with the 
environment causes the initial pure state of the central 
system to evolve into a mixed state, described by a re- 
duced density matrix (35j , obtained by tracing out all the 
degrees of freedom of the environment @, H, 0, EH • If the 
Hamiltonian of the central system H c is a perturbation, 
relative to the interaction Hamiltonian H ce , the pointer 
states are eigenstates of H ce [l5|, [3l[ . On the other hand, 
if H ce is much smaller than the typical energy differences 
in the central system, the pointer states are eigenstates 
of H c , that is, they may be singlet or triplet states. In 
fact, as we will show, the selection of the eigenstate as 
the pointer state is also determined by the state and the 
dynamics of the environment. 

In the simulations that we discuss in the paper, the 
interactions between the central system and the environ- 
ment are either Ising or Heisenberg-likc. The interest- 
ing regime for decoherence occurs when each coupling of 
the central system with the environment is weak, that 
is, A <C |J|, but there is of course nothing that pre- 
vents us from performing simulations outside this regime. 
The interaction within the environment are taken to be 
Heisenberg-like, f2 being a parameter that we change. 
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FIG. 1: (Color online) The time evolution of the correlation 
(>&(i)|Si ■ S2j 1 3'(i)) of the ferromagnetic central system with 
Heisenberg-like H ce and H a . The model parameters are A = 
0.15 and a: Q, = 0.075; b: fl = 0.15; c: Q = 0.20; d: fl = 0.30; 
e: SI = 1. The number of spins in the environment is N = 14. 




t\J\ 



FIG. 2: (Color online) The time evolution of the diagonal 
matrix elements of the reduced density matrix of the central 
system for A = 0.15 and Q. = 0.15 (case (b) of Fig. [[]). The 
number of spins in the environment is N = 14. 



Clearly, in case (a), the relaxation is rather slow and 
confirming that there is relaxation to the ground state 
requires a prohibitively long simulation. For cases (b) - 
(d), the results are in concert with the intuitive picture of 
relaxation due to decoherence: The correlation shows the 
relaxation from the up-down initial state of the central 
system to the fully polarized state in which the two spins 
point in the same direction. 

An important observation is that our data convincingly 
shows that it is not necessary to have a macroscopically 
large environment for decoherence to cause relaxation to 
the ground state: A spin-glass with N — 14 spins seems 
to be more than enough to mimic such an environment. 
This observation is essential for numerical simulations of 
relatively small systems to yield the correct qualitative 
behavior. 

Qualitative arguments for the high efficiency of the 
spin-glass bath were given in Ref. [22|. Since the spin- 
glasses possess a huge amount of the states that have an 
energy close to the ground state energy but have wave 
functions that are very different from the ground state, 
the orthogonality catastrophe, blocking the quantum in- 
terference in the central system [Til ITa ] is very strongly 
pronounced in this case. 

This conclusion is further supported by Fig. [2] where we 
show the diagonal elements of the reduced density matrix 
for case (b). After reaching the steady state, the nondi- 
agonal elements exhibit minute fluctuations about zero 
and are therefore not shown. From Fig. [21 it is then clear 
that central system relaxes to a mixture of the (spin-up, 
spin-up), (spin-down, spin-down), and triplet state, as 
expected of intuitive grounds. In case (e), the character- 
istic strength of the interactions between the spins in the 
environment is of the same order as the exchange cou- 
pling in the central system (f2 ss J), a regime in which 
there clearly is significant transfer of energy, back-and- 
forth, between the central system and the environment. 

^From the data for (b) - (d), shown in Fig.Q] we con- 
clude that the time required to let the central system re- 
lax to a state that is close to the ground state depends on 
the energy scale (fi) of the random interactions between 
the spins in the environment. As it is difficult to define 
the point in time at which central system has reached its 
stationary state, we have not made an attempt to char- 
acterize the dependence of a relaxation time on Q. 



III. HEISENBERG-LIKE H ce 
A. Ferromagnetic central system 

In this section, we consider a ferromagnetic (J = 1) 
central system that interacts with the environment via 
a Heisenberg-like interaction (recall that throughout this 
paper the environment itself is always Heisenberg-like) . 

In Fig. [TJ we present simulation results for the two-spin 
correlation function for different values of the parameter 
£1 that determines the maximum strength of the coupling 
between the N(N—l)/2 pairs of spins in the environment. 



B. Antiferromagnetic central system 

We now consider what happens if we replace the fer- 
romagnetic central system by an antiferromagnetic one. 

The main difference between the antiferromagnetic and 
the ferromagnetic central system is that the ground state 
of the former is maximally entangled (a singlet) whereas 
the latter is a fully polarized product state. 

In Fig. [31 we present simulation results for the two-spin 
correlation function for different values of the parameter 
f2. In passing, we mention that in our simulations, we 



5 



A 
-»— » 

CM 

C/J 



CO 

V 




-1—' 

CN 

W 



CO 

V 



i — 1 — i — ' — r 




J i I i L_ 



_l i I i L_ 



t|J| 



1000 2000 3000 4000 5000 6000 7000 8000 9000 
t|J| 



FIG. 3: (Color online) The time evolution of the correlation 
($(t)|Si ■ S2 1 of the antiferromagnetic central system 
with Heisenberg-like H ce and H e . The model parameters are 
A = 0.15 and a: fl = 0.075; b: Q = 0.15; c: Q = 0.20; d: 
Q, = 0.30; e: fi = 1. The number of spins in the environment 
is TV = 14. 



FIG. 5: (Color online) The time evolution of the correlation 
($(t)|Si ■ S2 1 of the antiferromagnetic central system 
with Ising-like H ce and Heisenberg-like H e . The model pa- 
rameters are A = 0.075 and a: fi = 0.075; b: Q = 0.15; c: 
Q = 0.30; d: fl = 1. The number of spins in the environment 
is TV = 16. 
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FIG. 4: (Color online) The time evolution of the diagonal 
matrix elements of the reduced density matrix of the central 
system for A = 0.15 and O = 0.15 (case (b) of Fig. [3]). The 
number of spins in the environment is TV = 14. 



change the sign of J only, that is we use the same pa- 
rameters for H ce and H e as in the corresponding simula- 
tions of the ferromagnetic case. Apart from the change is 
sign, the curves for all cases (a-e) in Fig. Q] and Fig. [3] are 
qualitatively similar. However, this is a little deceptive. 

As for the ferromagnetic central system, in case (a) , the 
relaxation is rather slow and confirming that there is re- 
laxation to the ground state requires a prohibitively long 
simulation. In case (e), we have f2 « | J\ and as already 
explained earlier, this case is not of immediate relevance 
to the question addressed in this paper. For cases (b) - 



(d), the results are in concert with the intuitive picture 
of relaxation due to decohcrcnce except that the central 
system does not seem to relax to its true ground state. In- 
deed, the two-spin correlation relaxes to a value of about 
0.65 - 0.70, which is much further away from the ground 
state value —3/4 than we would have expected on the 
basis of the results of the ferromagnetic central system. 
In the true ground state of the whole system, the value 
of the two-spin correlation in case (b) is —0.7232, hence 
significantly lower than than the typical values, reached 
after relaxation. On the one hand, it is clear (and to be 
expected) that the coupling to the environment changes 
the ground state of the central system, but on the other 
hand, our numerical calculations show that this change 
is too little to explain the apparent difference with the 
results obtained from the time-dependent solution. 

In Fig. 21 we plot the diagonal matrix elements of 
the density matrix (calculated in the basis for which the 
Hamiltonian of the central system is diagonal) for case 
(b). From this data and the fact that the nondiagonal 
elements are negligibly small (data not shown), we con- 
clude that the central system relaxes to a mixture of the 
singlet state and the (spin-up, spin-up) and (spin-down, 
spin-down) states, the former having much more weight 
(0.9 to 0.05) than the two latter states. Thus, at this 
point, we conclude that our results suggest that decoher- 
ence is less effective for letting a central system relax to 
its ground state if this ground state is entangled than if it 
is a product state. Remarkably, this conclusion changes 
drastically when we replace the Heisenberg-like H ce by 
an Ising-like H ce , as we demonstrate next. 
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IV. ISING-LIKE H ce 

In our simulation, the initial state of the central system 
is HI) and this state has total magnetization M = 0. For 
Ising H ce with Heisenberg-like H e coupling, the magneti- 
zation M of the central system commutes with the Hamil- 
tonian {T]) of the whole system. Therefore, the magneti- 
zation of the central system is conserved during the time 
evolution, and the central system will always stay in the 
subspace with M = 0. In this subspace, the ground 
state for antiferromagnetic central system is the singlet 
state \S) while for the ferromagnetic central system the 
ground state (in the M = subspace) is the entangled 
state \Tq). Thus, in the Ising-like H ce , starting from the 
initial state |||), the central system should relax to an 
entangled state, for both a ferro- or antiferromagnetic 
central system. 

If the initial state of the central system is |f J.), it can 
be proven (see Appendix) that 



(tf (t)|S! • S 2 |*(i)) F + (*(t)|S! • S 2 \^(t)) A = --, (4) 
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FIG. 6: (Color online) The time evolution of the correlation 
(>&(£) |Si ■ S2 1 \&(£)) of the antiferromagnetic central system 
with Ising-like H CE and Heisenberg-like H e - a: A = 0.0375 
and ft = 0.15; b: A = 0.075 and fi = 0.15; c: A = 0.075 and 
£1 — 0.3; d: A = 0.15 and fi = 0.3. The number of spins in 
the environment is N = 16. 



V. ROLE OF A 



where the subscript F and A refer to the ferro- antiferro- 
magnetic central system, respectively. Likewise, for the 
concurrence we find Cp(t) = C'a (t) and similar symme- 
try relations hold for the other quantities of interest. Of 
course, this symmetry is reflected in our numerical data 
also, hence we can limit ourselves to presenting data for 
the antiferromagnetic central system with Ising-like H ce 
and Heisenberg-like H e . 

In Fig. [5l we present simulation results for the two-spin 
correlation function for different values of the parameter 
Q. Notice that compared to Figs.[I]-[U we show data for 
a time interval that is three times larger. For the cases 
(b,c), the main difference between Fig. [3] and Fig. [5] is 
that for the latter and unlike for the former, the central 
system relaxes to a state that is very close to the ground 
state. Thus, we conclude that the presence of a conserved 
quantity (the magnetization of the central system) acts 
as a catalyzer for relaxing to the ground state. Although 
it is quite obvious that by restricting the time evolution 
of the system to the M = subspace, we can somehow 
force the system to relax to the entangled state, it is by 
no means obvious why the central system actually does 
relax to a state that is very close to the ground state. 

Intuitively, we would expect that the presence of a 
conserved quantity hinders the relaxation and indeed, 
that is what we observe in cases (a,b) where the relax- 
ation is much slower than in cases (a,b) of Fig. [1] or of 
Fig. [3] Notwithstanding this, in the presence of a con- 
served quantity, the central system relaxes to a state that 
is much closer to true ground state than it would relax 
to in the absence of this conserved quantity. 



Now, we study the effect of changing the strength A 
of the coupling between central system and the environ- 
ment. For a qualitative discussion of this aspect, it suf- 
fices to consider the case of Ising-like H ce , as we have 
seen that then, the central system most easily relaxes to 
its ground state. 

In Fig. [5J we present some representative simulation 
results for the two-spin correlation function for different 
values of the parameters A and fl. By simply compar- 
ing the time intervals of the plots for cases (a,b) and 
(c,d), it is immediately clear that the speed of relaxation 
changes drastically with A. For a "slow" environment 
(small enough fl) the effect is rather trivial, namely, the 
larger A the faster the relaxation. In the case (c) the sys- 
tem comes close to the triplet state in comparison with 
(d), probably, since the perturbation of the ground state 
of the central system is smaller. 

VI. SENSITIVITY OF THE RESULTS TO 
CHARACTERISTICS OF THE ENVIRONMENT 

Finally, we study the effect of small changes to the 
initial state of the environment and of the number of 
spins in the environment. 

For the spin glasses, the true ground state is rather 
hardly reachable and there are a lot of states with a very 
close energy but essentially different characteristics. To 
check how relevant it can be for our observations, we re- 
place the environment ground state by one of such states 
and study the time evolution of the central system as 
we did before. In Fig. \7\ we show typical results for 
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FIG. 7: (Color online) The time evolution of the correlation 
(>&(£) |Si ■ S2 1 of a ferromagnetic central system with 
Heisenberg-like H ce and Heisenberg-like H e with A = 0.15 
and f2 = 0.3. Initial state of the environment is solid line (a): 
ground state; dashed line (b): close to but not the same as 
the ground state. The number of spins in the environment is 
N = 14. 



a ferromagnetic central system with Heisenberg-like H ce 
and Heisenberg-like H e . In the initial state, the energy 
of the environment Eb — —2.247, which is a little bit 
higher than the ground-state energy of the environment 
E a = —2.321. The time evolution of the correlation func- 
tion of the two central spins for the cases (a) and (b) (see 
Fig. [7]) clearly demonstrates that in both cases, the cen- 
tral system evolves to the ground state, and that the 
dynamics of this evolution is also very similar. This con- 
firms that as long as the energy of the initial state of the 
environment is close to its ground state energy, the qual- 
itative features of the decoherence process remain the 
same. If, on the other hand, we prepare the environment 
in a random state (which, roughly speaking, corresponds 
to a very high temperature), the central system does not 
relax to its ground state but to a mixed state with a 
diagonal density matrix, as expected (see Fig. |5J| . 

Second, we study the effect of finite size of the environ- 
ment on the decoherence process. Some typical results 
for a ferromagnetic central system with Heisenberg-like 
H ce and Heisenberg-like H e with different numbers of the 
environment spins N are shown in Fig. [9] It looks rea- 
sonable to define the border between a mesoscopic and a 
macroscopic environment as a value of N for which the 
oscillations in the two-particle correlation are no longer 
well-defined. Thus, on the basis of the data displayed in 
Fig. [5] one can say that N w 11 is large enough for the 
spin-glass environment to mimic the macroscopic system. 
Needless to say, this statement is very qualitative but, in 
any case, the N dependence of the results shown in Fig. [9] 
demonstrate the effectiveness of the spinglass as a model 
environment to study decoherence processes with rather 
modest requirements to the environment size. 



VII. SUMMARY 

We have presented the results of simulations that ad- 
dress the question how a small quantum system evolves 
to its ground state when it is brought in contact to an 
environment consisting of quantum spins. Our system- 
atic study confirms the suggestion of Ref. [22| that the 
use of spin-glass thermal bath is indeed a very efficient 
way to simulate decoherence processes. Environments 
containing 14 - 16 spins are sufficiently large to induce a 
complete decay of the Rabi oscillations, this in sharp con- 
trast to environments that have a more simple structure, 
such as spin-chains or square lattices [22j |. 

In general, it turns out that the relaxation to the 
ground state is a more complicated process that one 
would naively expect, depending essentially on the ratio 
between parameters of the interaction and environment 
Hamiltonians. Two general conclusions are: (i) the cen- 
tral system more easily evolves to its ground state when 
the latter is less entangled (e.g., up-down state compared 
to the singlet) and (ii) constraints on the system such as 
existence of additional integrals of motion can make the 
evolution to the ground state more efficient. 

At the first sight, the latter statement looks a bit coun- 
terintuitive since it means that it may happen that a 
more regular system exhibits stronger relaxation than a 
chaotic one. The reason that is may happen is that the 
larger the dimensionality of available Hilbert space for 
the central system is, the more complicated the deco- 
herence process is due to appearance of the whole hier- 
archy of the decoherence times for different elements of 
the reduced density matrix. A manifestation of this phe- 
nomenon has been observed earlier [iff]: Under certain 
conditions, the same central system as studied here (four 
by four reduced density matrix) displays "quantum oscil- 
lations without quantum coherence" whereas for a single 
spin in magnetic field (two by two reduced density ma- 
trix) decoherence can, relatively easily, suppress the Rabi 
oscillations completely. 

We believe that these results can stimulate further 
development and clarification of the "decoherence pro- 
gram" [15l.l36j|. Assuming that the interaction with an en- 
vironment is weak enough, a hypothesis that the pointer 
states should be the eigenstates of the Hamiltonian of the 
central system was proposed |3lj . with the very ambi- 
tious aim to explain the basic phenomenon of "quantum 
jumps". In this paper, we demonstrate that, apart from 
just a strength of different interactions, also their symme- 
try and the amount of entanglement of the ground state 
of the central system may play an essential role. Among 
the cases which we consider in this paper, there are two 
situations where the standard decoherence scenario works 
as envisaged [3l| . If the ground state is not entangled (as 
in the case of the up-down state for the case of ferromag- 
netic interactions) or if the Hilbert space is restricted 
due to some conservation laws (as for the singlet ground 
state in the Ising-type interaction Hamiltonian), the cen- 
tral system clearly evolves to its ground state, supposed 
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FIG. 8: (Color online) The time evolution of the diagonal elements (left panel) and the real parts of the off-digonal elements 
(right panel) of the reduced density matrix in the antiferromagnetic central system, with Heisenberg-like H ce and Heisenberg- 
like H e (A = 0.15 and Q — 0.15). The initial state of the central two spins is the up-down state, and the environment is initially 
in a random state. The number of spins in the environment is TV = 14. 
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FIG. 9: (Color online) The time evolution of the correlation 
(i ) | Si ■ S2 1 of a ferromagnetic central system with 
Heisenberg-like H ce and Heisenberg-like H e with A = 0.15 
and f2 = 0.3. The number of spins in the environment is a: 
TV = 8; b: TV = 9; c: TV = 10; d: TV = 11; e: TV = 12. 



to be the pointer state according to Ref. [3l| . However, if 
the ground state of the central system is the fully entan- 
gled singlet state, and interaction Hamiltonian is generic, 
without symmetries, the system evolves to some mixture 
of the ground state and excited states. Of course, the 
data presented here are not sufficient to make strong, 
general statements about the character of the pointer 
states but we hope that, at least, our work will stimulate 
further research to establish the conditions under which 
the conjecture that the pointer states are the eigenstates 
of the central system hold. 



Appendix 

For the Hamiltonian Eq. Q, if A^ = A^ = 0, H ce 
is Ising-like and it is easy to prove that [M, H] = 0, 
implying that the magnetization of the central two spins 
is a conserved quantity. In our simulations, we take as 
the initial state of the central system the spin- up - spin- 
down state (| TO = (\S) + \T ))/V2). Hence, because 
[M, H] = 0, the central spin system will always stay in 
the subspace of M = 0. Thus, at any time t, the state of 
the whole system can be written as 



m)) = \S)\</>s(t)) + \T )\<t>T {t)), (5) 

where \<j>g) and \4>t } denote the states of the environ- 
ment. 

Let us denote by {|V>i)} a complete set of states of 
the environment. Within the subspace spanned by the 
states {IS 1 ) IV'i)) l^o)!^)}) the Hamiltonian Eq. |T]) can be 
written as 

H = E s \S)(S\+E T \T )(T \+H e 
1 N 

-o E( A S - A B) (IWol + |T )(5|) II (6) 

where we used {S\Sl\S) = (T \Sf\T ) = {S\S$\S) = 
(T \SI\T ) = 0, (T \Sf\S) - 1/2, and (T |Sf |S) = -1/2. 

Introducing a pseudo-spin a = (p x , a v , <r z ) such that 
the eigenvalues +1 and — 1 of <j z correspond to the states 
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\S) and | T ), respectively, Eq. ([6]) can be written as 

E s - E T E s + E T 
H = + a +H e 

1 N 

-i]T(Ag-Ag)/;a*, (7) 
j=i 

showing that in the case of Ising-like H ce , the model 
Eq. ((T|) with two central spins is equivalent to the model 
Eq. J7]) with one central spin. 

^From Eq. ((7]), it follows immediately that the Hamil- 
tonian is invariant under the transformation {J, a z } — > 
{ — J, — <j z }. Indeed, the first, constant term in Eq. ((7| is 
irrelevant and we can change the sign of the second term 
by rotating the speudo-spin by 180 degrees about the x- 



axis. Therefore, if the initial state is invariant under this 
transformation also, the time-dependent physical proper- 
ties will not depend on the choice of the sign of J, hence 
the ferro- and antiferromagnetic system will behave in 
exactly the same manner. 

For the case at hand, the initial state can be written as 
( | S) + 1 To ) ) 1 4>o ) / V% , which is trivially invariant under the 
transformation a z — > — a z . Summarizing: For Ising-like 
H ce (A-^ — A- j — 0), and an initial state that is in- 
variant for the transformation IS") <-> |T )), (*(t)| A\^(t)) 
does not depend of the sign of J, for any observable A 
of the central system that is invariant for this transfor- 
mation. Under these conditions, it is easy to prove that 
Eq. (0| holds and that the concurrence does not depend 
of the sign of J. 
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